The Aquarius Salinity Retrieval Algorithm 
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Abstract — This paper gives an overview of the algorithm for re- 
trieving sea surface salinity from the AQUARIUS L-band radi- 
ometer and its physical background. 
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I. Background 

The Aquarius L-band radiometer/scatterometer system is 
designed to provide monthly salinity maps at 150 km spatial 
scale to an accuracy of 0.2 psu [1]. The sensor was launched 
on June 10, 2011, aboard the Argentine CONAE SAC-D 
spacecraft. The L-band radiometers and the scatterometer have 
been taking science data observations since August 25, 2011. 

II. Overview of the Aquarius Level 2 Processing 

This section gives an overview over the Aquarius salinity 
retrieval algorithm. 

The instrument calibration [2] converts Aquarius radiome- 
ter counts into antenna temperatures (TA). The salinity retriev- 
al algorithm converts those TA into brightness temperatures 
(TB) at a flat ocean surface and sea surface salinity (SSS). As 
a first step, contributions arising from the intrusion of solar, 
lunar and galactic radiation are subtracted. The antenna pattern 
correction (APC) removes the effects of cross-polarization con- 
tamination and spillover. The Aquarius radiometer measures 
the 3rd Stokes parameter in addition to vertical (v) and horizon- 
tal (h) polarizations, which allows for an easy removal of io- 
nospheric Faraday rotation. The atmospheric absorption at L- 
band is almost entirely due to molecular oxygen, which can be 
calculated based on auxiliary input fields from numerical 
weather prediction models (e.g. NCEP) and then successively 
removed from the TB. The final step in the TA to TB conver- 
sion is the correction for the roughness of the sea surface due to 
wind, which is addressed in more detail in section III. 

The TB of the flat ocean surface can now be matched to a 
salinity value using a surface emission model that is based on a 
model for the dielectric constant of sea water [3], [4] and an 
auxiliary field for the sea surface temperature. In the current 
processing only v-pol TB are used for this last step. 


III. Model for the Emission of Wind Roughened 
Ocean Surface at L-Band 

One of the major error sources in the Aquarius salinity re- 
trieval is the change in the ocean surface emissivity due to 
roughening by wind. 

A. Wind Speed Dependence 



Figure 1: Harmonic coefficients of the expansion (1) of the wind induced 
surface brightness temperature AE W • T s (in Kelvin) for the 3 Aquarius radi- 
ometers: blue = v-pol, red=h-pol. The Earth incidence angles for the beams 
are: 28.7° (RAD 1), 37.9° (RAD 2) and 45.5° (RAD 3). 

A model for the wind induced emissivity has been devel- 
oped post-launch based on actual Aquarius observations. As a 
first step, Aquarius brightness surface brightness temperatures 
were collocated with surface wind speed and wind direction 
measurements from NCEP, SSMIS F17 and WindSat, from 


which the harmonic coefficients A i? i=0, 1 ,2 of the wind in- 
duced emissivity AE W can be de derived: 

AE w (W,cp r )-T s = 

A 0 (w)+A>W)-Aos^ r )cos 2 ^ )• ( r ) 

W is the surface wind speed, cp r is the surface wind direction 
relative to the looking azimuth and T s is the sea surface tem- 
perature. The results for the 3 Aquarius horns are shown in 
Figure 1. This allows performing a roughness correction 
based on auxiliary wind speed and wind direction fields from 
NCEP. 



Figure 2: Comparison of the wind induced surface brightness temperature 
AE W • T s [• 290 Kelvin] from Aquarius (red) with the results of: [5] (black), 
[6] (green) and [7] (blue). 


Figure 2 compares the L-band wind induced surface bright- 
ness temperature as function of wind speed from Aquarius 
with other measurements. 

B. Roughness and Scatterometer Backscatter 



Figure 3: The wind induced surface brightness temperature AE W • T s [• 290 
Kelvin] binned with respect to surface wind speed from NCEP W NC ep and the 
scatterometer measurement o 0 ,w- The o 0 ,w has been scaled to units of wind 
speed [m/s]. The wind direction signal has been removed from both AE W • T s 
and Oo,w • 


As a next step the radar backscatter cross section measure- 
ments a 0 of the Aquarius scatterometer can be used in addition 
to NCEP wind speed and direction in order to improve the 
roughness correction. Both, radiometer and scatterometer 
observation need to be first corrected for wind direction, 
which is done based on the functions from Figure 1 and an 
analogous analysis for o 0 . The wind induced emissivity can be 
parameterized and calculated from a lookup table containing 


both NCEP wind speed W NC ep and the g 0 for the vv-pol back- 
scatter. The results are shown in Figure 3. 


C. Surface Roughness Correction 



Figure 4: Schematic flow of the surface roughness correction algorithm. 


Figure 4 shows the flow of the algorithm that correcting for 
the surface roughness. If the [W NCEP , o 0 ,vv] bin in the 2- 
dimensional lookup table (Figure 3) is underpopulated, only 
W NC ep is used in the roughness correction. 

IV. First Results for Sea Surface Salinity 



Figure 5: Aquarius salinity map for the month of September 2011. 



Figure 6: Statistics of Aquarius versus HYCOM salinity (SSS) as function of 
NCEP wind speed: dashed line = bias, full lines = standard deviation. 
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Figure 7: Statistics of Aquarius versus HYCOM salinity (SSS) as function of 
SST: dashed line = bias, foil lines = standard deviation. 
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